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Abstract: One-loop quantities in QFT can be computed in an efficient way us- 
ing the worldline formalism. The latter rests on the ability of calculating ID path 
integrals on the circle. In this paper we give a systematic discussion for treating 
zero modes on the circle of ID path integrals for both bosonic and supersymmetric 
nonlinear sigma models, following an approach originally introduced by Friedan. We 
use BRST techniques and place a special emphasis on the issue of reparametrization 
invariance. Various examples are extensively analyzed to verify and test the gen- 
eral set-up. In particular, we explicitly check that the chiral anomaly, which can be 
obtained by the semiclassical approximation of a supersymmetric ID path integral, 
does not receive higher order worldline contributions, as implied by supersymmetry. 
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1. Introduction 

The worldline formalism is an efficient and economical way of calculating Feynman 
diagrams |[l|]. It describes the propagation of various relativistic particles by one 
dimensional (ID) path integrals. Recently we have shown how to extend this method 
to the case of spin and spin 1/2 particles coupled to background gravity 0, Q . Other 
applications run from the calculation of the heat kernel to the calculation of chiral 
and trace anomalies. 

In all these cases, the ID path integrals are calculated on the circle with a fi- 
nite propagation time, the proper time. In many circumstances this proper time is 
integrated over, as it represents the only modulus of the circle. The actions in the 
path integrals are those of one dimensional nonlinear sigma models. They describe 
the propagation of particles in a curved background, just like 2D nonlinear sigma 
models describe the propagation of strings^. Sigma models are super-renormalizable 
in ID, and their UV structure together with the necessary renormalization condi- 
tions (which produce explicit counterterms) have been extensively discussed in the 
literature using various regularization methods 0, |[, ||. 

^The requirement of conformal invariance restricts the possible backgrounds on which the string 
propagates, but no such requirement is present for the particle case, at least for spin and 1/2. 
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In this paper we plan to address in a more systematic way the infrared issue 
related to the treatment of zero modes appearing on the circle and the interplay 
with the reparametrization invariance of nonlinear sigma models. 

Different ways of treating such zero modes have been developed in the literature. 
For example in ^ an arbitrary background charge function p(r) was used to inter- 
polate between various boundary conditions, thus lifting the zero modes in different 
ways. There it was shown that for models in flat space the effective lagrangian cal- 
culated with different background charges p differed only by total derivatives, thus 
producing the same effective action. While this is not causing any particular problem 
in flat space (total derivatives are even beneficial in certain cases, since they allow 
to cast the effective action in a more compact form), the naive extension of this 
method to curved space was seen to introduce noncovariant total derivatives p, 0], 
raising questions about the correct use of Riemann normal coordinates to simplify 
calculations. 

A general method for dealing with these zero modes has been developed by 
Friedan in his treatment of 2D bosonic sigma models and recently employed in 
the ID case by Kleinert and Chervyakov |T^, where the comparison between the so- 
called SI (string inspired) propagator and the DBC (Dirichlet boundary conditions) 
propagator was carried out to show how the former could produce a covariant result 
as the latter. The string inspired propagator |]lT| is translational invariant on the 
worldline, while the DBC one is not. 

Here we review and analyze the treatment of the zero modes to clarify it further 
and resolve some remaining puzzles. We use ERST methods to factor out the zero 
modes. The BRST symmetry is related to the gauge fixing of a shift symmetry. It is 
the shift symmetry typical of background field methods [^. However, in the present 
case the "background field" is integrated over in the path integral (i.e. it is made 
dynamical), and thus the shift symmetry must be gauge fixed. An immediate result of 
this procedure of extracting the zero modes is that different gauges are guaranteed 
to produce the same effective action, implying that effective lagrangians can only 
differ by total derivatives, even in curved space. However, the total derivatives 
will in general be noncovariant. On the other hand, suitable covariant gauges will 
naturally produce effective lagrangians differing from each other only by covariant 
total derivatives. Various examples are extensively analyzed to test these predictions. 
In particular, we show how the calculation of the trace of the heat kernel and the 
related Seeley-DeWitt coefficients, which identify the effective action for a scalar 
particle in the proper time expansion [0, is achieved with different background 
charges. 

We also extend the zero mode treatment to the supersymmetric case. This is 
relevant for the worldline description of spin 1/2 fermions coupled to gravity . We 
use it for an explicit check that the chiral anomaly does not receive higher order 
worldline contributions. In fact, one may recall that the = 1 supersymmetric 
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nonlinear sigma model was used in [14| to compute the chiral anomaly of a spin 1/2 
field. The computation was based on the fact that the chiral anomaly could be iden- 
tified as the Witten index of the corresponding supersymmetric quantum mechanical 
model [15|. Supersymmetry implies that higher order worldline contributions should 
not modify the value of the Witten index, and we test this explicitly using ID path 
integrals. 

The paper is organized as follows. In section 2 we discuss bosonic nonlinear sigma 
models. As a test we calculate perturbatively the trace of the heat kernel and the 
related Seeley-DeWitt coefficients with different background charges. In section 3 
we consider supersymmetric nonlinear sigma models and use them to study worldline 
corrections to the susy quantum mechanical computation of the chiral anomaly. In 
section 4 we present our conclusions. For completeness and further clarifications we 
present the simpler case of a bosonic linear sigma model in appendix |A.1| . Other 
appendices include conventions and integrals needed for the computations described 
in the text. 



2. Bosonic nonlinear sigma models 

Let us consider the partition function of the ID nonlinear sigma model 

(2.1) 



Z(/3) = (pVx e-^["] , S[x] = ^ j^dr 



where Qf^u and V are a metric and a scalar potential defined on target space, which 
we take to be D dimensional. The path integral is computed with periodic boundary 
conditions (PBC), i.e. on the circle. We use euclidean time, and it is well-known 
that periodic boundary conditions then yield the statistical partition function. The 
circle is just the loop made by the particle in target space. It can be parametrized 
by t G [0, P], with /3 the total length of the circle. We also use a rescaled proper time 
T = t/P and this rescaling explains the factor ^ multiplying the action as well as the 
factor in front of the scalar potential. 

The partition function Z{(3) is sometimes called the trace of the heat kernel^ [0 



It is related to the one-loop effective action of a scalar field with kinetic operator 
— □ + 2V + ( □ is the covariant scalar laplacian depending on the metric g^u and 
m is the mass of the scalar particle) by an integral over the proper time (3 

n9,V] = -lJ^^^e-'^-'^Z{f3). (2.2) 

In a perturbative computation zero modes appear. In fact for g^^,^ = Sfj,^ and 

V = the action is invariant under the constant translations (5x'^(r) = e'*. Hence 

^The operator e~^^ is called the heat kernel, and the partition function is given by Z{f]) — 
Tr e^^^ = <f Vx e~^^^\ where iS'[a;] is the action corresponding to the model with quantum hamil- 
tonian H. 
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the volume of target space appears as a factor, just like the volume of a gauge 
group. Background fields generically break this translation invariance. Nevertheless 
it is both useful and necessary to extract from the path integral the "collective 
coordinates" or "center of mass" of the loops a;^(r), which for simplicity we continue 
to call zero modes. It is useful, since it allows to produce the partition function 
as an integral of a partition function density. It is necessary, since in perturbative 
calculations around the free action one needs to invert the free kinetic term to obtain 
the perturbative propagator. A general method for treating the zero modes for 
nonlinear sigma models has been developed by Friedan [Q, and employed recently in 
the ID case by Kleinert and Chervyakov [|lOl. Useful references are also |jl6|, |17| . 

Let us rederive these results and extend them by using an arbitrary background 
charge p{t). To extract the zero modes one can proceed as follows^. The action of the 
nonlinear sigma model S'[a;(r)] depends on the periodic paths x^{t) which describe 
loops with the topology of a circle in target space. One may introduce a redundant 
variable Xq by setting x^{t) = Xq + ?/^(r) in the action, S'[x(r)] = S[xo + ?/(t)]. Of 
course this automatically introduces the shift symmetry 







e 



6y^{T) = -e^ (2.3) 



as in the background field method [|T2|. We only consider constant Xq so that the shift 
symmetry requires a constant parameter e^. However, contrary to the background 
field method, we now consider both Xq and y^ir) as dynamical variables (i.e. to 
be integrated over in the path-integral). The shift symmetry is thus promoted to a 
gauge symmetry and must be gauge fixed since each physical configuration has to 
be counted only once. To gauge fix we use BRST methods and introduce a constant 
ghost field r]'^ together with the following BRST transformation rules 

6x^ = r^^A 

67]^" = . (2.4) 

To fix a gauge one must also introduce constant nonminimal fields fj^, tt^ with the 
BRST rules 

= in^A , (57r^ = . (2.5) 

Choosing the gauge fixing fermion 

^ = % [ dTp{T)y^^{T) (2.6) 

^0 



3t 



In appendix A.l we describe the simpler case of a linear sigma model using both Faddeev-Popov 



and BRST methods. 
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which depends on the arbitrary function p(r) normahzed to drpij) = 1 produces 
the following gauge fixed action 



Sgf[xo,y,r],r],n] = S[xo + y] + — ^ 




where ^ denotes a BRST variation with the anticommuting parameter A removed 
form the left. In this gauge the ghosts can be trivially integrated out, while the 
integration over the auxiliary variable produces a delta function which constrains 
the fields y^ to satisfy 



With this constraint the perturbative kinetic term for the periodic fields y'^ir), pro- 
portional to can be inverted to obtain the propagator. The BRST symmetry 
implies that the partition function is independent of the gauge parameter p. The 
specific case of p(r) = 6{t) gives the DEC propagator since y^{0) = y'^{l) = 0. The 
case p(r) = 1 gives instead the SI propagator since now the center of mass is absent 
from the fluctuations y^, see e.g. the discussion in |^. 

Thus the partition function is independent of p and can be expressed as an 
integral over the zero modes 



where the factor ^'L has been extracted for convenience from the definition of 



z^''\xo,P). Note however that the density z^'^\xo,P) may in general depend on p. 
This can only happen through total derivatives which must then integrate to zero. 

Although this is correct, at least formally, there are some practical problems. 
The constraint arising from the gauge fixing does not have simple transformations 
rules under change of coordinates (coordinate differences like = x^ — Xq do not 
transform as vectors). This causes some technical problems when one wants to 
check the explicit covariance of the final result. In particular, one explicitly finds in 
, 13) total derivatives which depend on the choice of p and are not covariant 
under change of the coordinates Xq H, 0] . Let us recall that DEC are related to the 
calculation of the heat kernel even for non coinciding points, and they are known 
to give a covariant result. The other propagators related to different background 
charges produce instead noncovariant total derivatives. A calculation at order f3 
using both the DEC and SI propagators was presented in to exphcitly identify 
the noncovariant total derivative term appearing at that order. For general p the 




(2.8) 




(2.9) 
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expression given in generalizes to 



ixo,P) = 1 + (l^^ - ^) + ^(^p + ^)^.(v^^?"' V) + OiP') (2.10) 



where the precise value of Cp is defined later in eq. (p.29|) . The appearance of 



noncovariant total derivatives may raise doubts about the use of Riemann normal 
coordinates which are often used to simplify calculations. In fact, in ||^ the assump- 
tion of a naive use of Riemann normal coordinates was seen to produce the wrong 
trace anomaly in 2 and 4 dimensions. 

Riemann normal coordinates can nevertheless be used, as showed by Friedan in 
his discussion of nonlinear sigma models 0. Friedan noticed that the simple linear 
shift symmetry ( |2.3| ) becomes nonlinear when using Riemann normal coordinates C,^ 
centered at Xq. One should use this nonlinear shift symmetry to correctly perform 
the gauge fixing in Riemann coordinates. Riemann normal coordinates have the 
property that they are manifestly covariant under reparametrization of the point Xq. 
Now the action S[xo,^{t)] = S[xo + y{xo,C,{T))] is invariant under the nonlinear shift 
symmetry 



Sxi^ 







e 



Se{r) = ~Q'uixo,aT))e^ (2.11) 



which is a reformulation of ( p. 3D in these new coordinates. Note that since the origin 
of the Riemann normal coordinates is shifted, the nonlinear transformation is defined 
in such a way that the new fields C,^' = + 6^^ are expressed in Riemann normal 
coordinates defined around the new origin Xq' = Xq+ 6xq. The expression Q'^uixo,^) 
can be explicitly calculated and can be found in p. We report it here up to the 
order needed in subsequent calculations 

+ (lv,V5i?^/5. - ^R''apxR\s.)C^^C^' + Oie) ■ (2.12) 

Now we can introduce ghosts and auxiliary fields as usual. The BRST symmetry 
for the nonlinear shift symmetry is 



6x^ 


= r]''A 




= -Q^ 




= 




= iTTfiA 




= . 



(2.13) 
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It is nilpotent since Q'^„{xo,^{t)) satisfies certain relations arising from tlie abelian 
nature of tlie sliift symmetry. Using tlie gauge fermion 

^ = rj, ['dTp{T)e{T) (2.14) 
Jo 



produces the gauge fixed action 



Sgf[xo, V, V, tt] = S[xo, ^] + — ^ 



Jo 

-ri^ [ dTp{T)Q>',ixo,ar))v'' ■ (2-15) 
Jo 



The integration over tfie auxiliary variable gives again a delta function which 
constrains the fields to satisfy 

fdTp{T)i'^{T)^Q (2.16) 

Jo 

so that their propagator can be obtained. This constraint has a simple tensorial 
transformation law under the change of coordinates of Xq, in fact these coordinates 
transforms as vectors under a reparametrization of the origin Xq. The ghosts now 
give a nontrivial contribution, i.e. a nontrivial Faddeev-Popov determinant. 

The above gauge fixed actions can be used in the path integral. Of course, one 
also needs to use a path integral measure that is both reparametrization and BRST 
invariant. This is given by 



Vx = Y\ V g{x{T))dx{T) — > dxodr]dr]d7r Y\ V 9{x{T))dx{T) 

0<T<1 0<T<1 

^ dxo df) dfj diT Y\ Vgi^o + y{T))dy{T) . (2-17) 

0<r<l 

where we have first added to the sigma model measure the measure for the BRST 
quartet x^.Tj^^fj^.Ti^ which is formally identical to unity (two commuting fields give 
a volume which is compensated by that of the two anticommuting fields), and then 
performed the change of variables identified by the "background-quantum" split x^ — 
Xq + y^. We have here used the linear splitting, but we could equally well use the 
nonlinear one written in terms of the Riemann normal coordinates centered at Xq. 
In fact the measure is reparametrization invariant, i.e. of the same form in any 
coordinate system 

Vx ^ dxodrjdfjdTT \^9{x^^I(j))d^{r) . (2.18) 

0<T<1 
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Note that with g{xo,^{T)) we have indicated the determinant of the metric in Rie- 
mann normal coordinates centered at xq. For future reference we hst the expansion 
of the metric in Riemann coordinates centered at Xq up to the order needed in later 
calculations 

+ (^V5V,i?^„;3.(xo) + ^i?^./i?.^5.(a;o))reW + 0(e^) . (2.19) 

It is now useful to introduce additional ghosts to exponentiate the nontrivial 
part of the measure. We use commuting and anticommuting 6^, c'^ ghost fields to 
reproduce the correct measure 

n V9{^o,arm{r) = n ^^(^) / n da{T)db{r)dc{r) e-'-- (2.20) 

0<T<1 0<T<1 0<T<1 

S^sr[La,b,c\ = ], f dr Hg.Axo^Oia'ci'' + b^C)] . (2.21) 
The extra vertices arising from the measure will contribute together with similar 



vertices from the sigma model action to make the final result finite |19 



We are now ready to re-assemble all parts of the path integral with the zero 
modes factored out by the nonlinear shift symmetry 

Z{p) = j dxodr]dr]dn ^ D^DaL'6Dce-^''^ [^0'«'^'^'^l-^'"=^[^0'«''''''''=l . (2.22) 

The auxiliary field can be integrated out to obtain 



where 



j dxo dr] df] j DiDaDhDc 5(^j^ dr p{t)^^'{t)^ e^^" (2.23) 



= Sgf\xo, ^, r], rf\ + S„,sr[xo, a, b, c] (2.24) 

with the auxiliary field tt^ eliminated from 5*^/. Finally, when perturbation around 
the leading terms in (|2.12|) and (p.l9|) is appropriate, one immediately obtains the 
following perturbative expansion 



Z{(3) = I d^x,^^ (exp(-5(-))) (2.25) 

where the expectation value of the interactions are computed with the propagators 

{e{r)C{<y)) = -Pg'''{x,)B^,){T,a) 
{a^^{T)a''{a)) = (3g^'{x,)/\gH{r - a) 
{b^{ry{a)) = -2(3g^''{xo)A,,{r-a) 

iv'V.) = ■ (2.26) 
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The terms in the measure can be traced back as follows: the factor \/ g{xo) is due to 
the a, b, c ghosts which contain the constant "zero modes", the factor {2iTf3)~~ is the 
usual free particle measure which corresponds to the determinant oi —^-^ on the 
circle with zero modes excluded. The Green functions appearing in the propagators 
are as follows. i3(p) (r, a) is the Green function of the operator ^ acting on fields 
constrained by the equation dTp{T)^^{T) = 0. It depends on p and satisfies 

|^i3(,)(r,a)=5(r-a)-p(r) . (2.27) 
It is explicitly given by [0] 

i3(p)(r, or) = A(r - a) - Fp(r) - F,((t) + Cp (2.28) 



where 



11 1 

A(r — a) = -|r — crl (r — cr)^ 

^ ^2' '2^ ^12 

F,(r) = [ dx A(r - x)p{x) , = [ dr Fp{t)p{t) , (2.29) 



and it clearly satisfies 



/ dr p{t) i5(p)(r,o-) = . 
Jo 



(2.30) 



Note that the auxiliary function A(r — a) is the unique translational invariant Green 



function on the circle (the "string inspired" propagator of [|TT|, which corresponds to 
p(r) = 1). In the following we will simply write B(^p^ = as no confusion can arise. 
The Green function for the ghosts is given by the delta function 

Aghir -a)= Sir - a) (2.31) 

but we continue to call it Agh as in perturbative calculations it appears in a regulated 
form. 

We are now ready to test this set up. To summarize, the expectations are that the 
partition function Z{f3) will not depend on p, but that the density (5) will in 

general be p-dependent through total derivatives which integrate to zero. Moreover, 
using Riemann normal coordinates and the associated nonlinear shift invariance to 
extract the integral over the zero modes Xq, one expects z'^p\xq, 13) to be covariant 
under change of coordinates of Xq, so that the specific total derivatives which may 
eventually appear will also be covariant. 

In the next subsections we explicitly verify, using Riemann normal coordinates, 
that these total derivatives term are non-zero and covariant also in the presence of 
external potentials like V . 
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2.1 Partition function at 3 loops 

We present here the exphcit perturbative calculation of the partition function den- 
sity to order (3"^ using Riemann normal coordinates (RNC) and dimensional regu- 
larization on the worldline. The nonlinear sigma model in one dimension is super- 
renormalizable and one needs to choose a specific regularization scheme to compute 
unambiguously the perturbative expansion. We use dimensional regularization which 
requires an explicit counterterm Vdr = —\R to guarantee that the sigma model in 
will have H = —^D+V as quantum hamiltonian 0. In the following we will use 
the rules for dimensional regularization explained in Dimensional regularization 
has also been discussed for ID nonlinear sigma model with infinite proper time in 



The partition function density in eq. ( p.9| ) and ( p.25| ) can be expressed in terms 
of connected worldline graphs as 

/'\xo,P) = (exp(-5f = exp((e-^'"^ - 1) (2.32) 

where Sg = Sgj + Smsr is the full quantum action, {...)c denotes connected graphs, 
and the propagators are given in eq. (p.26|) . 



In order to appreciate the contribution from the FP determinant, we separate 
the corresponding action 

Sfp = -% [ dr p{r)Q':{xo, ^(r)) (2.33) 
Jo 

from the other contributions, gathered in S = Sg — Spp. As usual, we organize the 
interaction terms in such a way that Sg^^^ = Sg — Sg^2 = Sg^^ + Sg^^ + . . ., where 
Sg,n = Sn + Spp^n givcs risc to vevs of order P^~^ (recall that Sg^^ = in RNC, see 
eqs. ( |2A^ ) and {^J^ ). Hence 



)c - 1 = ( - 5,,4 - Sg,, + -Si,} + O{0') . (2.34) 



2 

We denote by Vg the potential which includes the counterterm arising in dimensional 
regularization 

V<i = y + Vdr = V -\r. (2.35) 



We now list the results for the various terms appearing in (|2.34| ) and report 



in the appendix |A.3| the expressions and values of the relevant connected one- and 



two- loop worldline integrals ifj, calculated in dimensional regularization 
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2 ('^FP,4)c 



{Sfp,^ 5*4) c 



1 

20 



-H^ + H,)nR+-{H,-H,) 
45 



p2 I _ d2 



+ 2^1 °K 



2C 



" 6 144; 



20 45 '^'^ 



I _ p2 

30 '^'^"^ 



+ 



1 

12 



□K 



36 



+ -(i7l3 + i^l4 



+ 

+ 
+ 



Cn H En 

162 ^ 27 



7 



12960 
13 



F?2 



HiQ dR 



30 



Hip. H -f^l7 

30 ^ 12 



45 ^^18 ^ 



d2 



10^" 



180 15 ^ 



540 " 



45 ^ 



4320 
7 



- —Cp + -D, 
60 ^ 180 ^ 30 ^ 



12960 
1 



R 



1 



9 



—Hij R^^ — — 



ci 



12 ^ 2880 



1 



R2 



2D, + E, + —] R 



9 
9 



18 



19 



2Hoa) R 



-2Cl 



3D. 



2 ^ 



720 



1 



180 



(2.36) 



All tensors appearing here are evaluated at Xq. 

As a check on these results, note that all contributions from the FP determinant 
vanish when using DBC. In fact, setting p(r) = b(r) one obtains for the p dependent 
coefficients C^, Dp, Bp the following values Cdbc = - 1^, Ddbc = Edbc = -lio- 
This is expected, since for p(r) = 5{t) the constraint (|2.16|) enforces the DBC, i.e. 
.^^(0) = ^'^(l) = 0. This in turn implies that the FP determinant in eq. ( p.l5| ) 
becomes trivial, since Q^u{,xq,Q) = 6^, as seen form eq. ( p. 121) . Note also that we 
have defined the p-dependent coefficients to vanish in the SI case, i.e. when p(r) = 1. 
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We can now insert the resuhs ( p.36 ) into ( 2.34 ) and ( 2.32| ) to obtain the partition 
function density vahd to the order 



Cr, — Dn H En 

144 ^ 12 24 
J20 720 ^"^""^ 



864 



dR 



(2.37) 



Expanding the exponent to order one sees that only covariant total derivatives 
have p-dependent coefficients 



z^''\xo,l3) = 1-/3 











+ 13^ 






2 V24 7 



_|_d2 3— R2 

720 ^"""^ 720 




1728 



i(c,-i)nv; 



(2.38) 



Let us discuss some consequences of this formula. We can integrate this density 
over the zero modes xq to obtain the partition function 



d^Xo 



(27r/3)f 



1 + ai{xo)f3 + a2(xo)/3^ + ... 



(2.39) 



where a„(xo) are the so-called Seeley-DeWitt coefficients. Inserted into ( ^.2| ) this 
gives in turn the effective action of the scalar field with kinetic term —n+2V+m'^. It is 
well-known that this proper time expansion of the effective action fails for massless 
fields. In fact for vanishing mass the damping factor e~2™'^'^, which guarantees 
convergence in eq. ( p.2|) , becomes unity. Nevertheless the Seeley-DeWitt coefficients 
are still useful in this well, as they give the counterterms needed to renormalize 

the one-loop effective action. Moreover, for conformal fields in D dimensions the 
coefficient ao gives the local trace anomaly |jT3[ . 
The main points to stress are: 

• For the validity of the perturbative calculation of Z{f3) we have to assume 
that all external fields describing the interactions should vanish sufficiently fast at 
infinity. Thus all total derivatives integrate to zero. Therefore gauge independence 
is verified. 

• The covariant local expansions of z^^^xq, P) in Riemann normal coordinates 
(RNC) are different for different p's. The difference is given by total derivatives with 
coefficients depending on p and which are explicitly nonvanishing. From the compu- 



tation to order P in this could not be evinced. In fact at that order there is no 
covariant total derivative with the correct dimensions that could possibly contribute. 
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Thus, at order (3 different p's produce the same local expansion. In principle there 
could have existed hidden relations guaranteeing the same local expansion of z{xo, /?) 
for different p's at any order in /3. We see explicitly that this is not the case and 
covariant total derivative may arise. 

• One can use these results to compute trace anomalies, which in D dimensions 
are given by the Seeley-DeWitt coefficient ao associated to the corresponding con- 
formal operator. In D = 2 there cannot be any covariant total derivative contribution 
to ai, so the local trace anomaly can be obtained with any p. In D = 4 there is a 
difference. However, it only affects the nR term, which is a trivial anomaly (it can 
be canceled by a counterterm) . Thus also in this case any p is good enough for the 
computation of the trace anomaly. To be specific, a. D = 4 conformal scalar needs 
a potential V = where ^ = ^^"^-^^^ = |. Therefore Vg, which includes the DR 
counterterm, becomes Vg = —^R- Inserting Vg into ( p.38| ) produces the following 
trace anomaly for a conformal scalar field 

1„ 1_ 1\„ lo 1 



which is the correct value as far as the universal terms are concerned. In fact, only 
the total derivative term has a p-dependent coefficient; the well-known value ^ is 
recovered with DEC 

• In it was argued that one could reach covariant results for the partition 
function density also using the "string inspired" propagator in arbitrary coordinates, 
i.e. through a noncovariant expansion. However the method proposed in section 
7 of |Ty] does not seem to be correct, as far as we understand it. In particular, if 



one were to include an external potential V in the path integral, there would be no 
terms in that method that could covariantize it. In any case, these facts may leave 
some doubts about the correct covariantization obtained using RNC in the presence 
of external potentials. To make sure that the previously described gauge-fixing in 
RNC achieves the correct covariantization even with external potentials, in the next 
subsection we present a four-loop computation in RNC which tests this issue. 

2.2 A look at order {3^ 

In the last subsection it was pointed out that RNC, with the correct factorization 
of the zero modes, should produce automatically the correct covariant result for 
the partition function density ,(3), also in the presence of terms that are not 

explicitly constructed from the metric tensor such as an external scalar potential V. 
In fact, gauge independence of Z{(3) guarantees that j3) should always contain 

a universal (p-independent) covariant part and a p-dependent total derivative, which 
vanishes upon space-time integration. The p-dependent total derivative is in general 
a noncovariant expression (as seen in eq. (|2.10| )) if calculated in arbitrary coordinates. 



It is instead covariant if calculated with RNC. 
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To make sure that the correct covariantization happens in RNC also in the 
presence of an external potential V, we compute the terms linear in V that arise 
at order They belong to the Seeley-DeWitt coefficient 03, and are enough to 
test the correct covariantization of the result. This test is in the same spirit of 
which employed RNC to a higher perturbative order to discover problems in the 
computation of the trace anomaly with the SI propagator. 

Thus we analyze at order f3^ the partition function density 



-.-S, 



(int) 



(2.41) 



where, as before, we have included in Sg*"*^ an arbitrary external scalar potential 
Vq = V — |-R. It is enough to consider the terms linear in Vq. As a consequence, we 
only need the perturbative expansion of the ghost action ( p.33| ) to order We thus 
only need to calculate the correlators 



-1 pI 

pVaVq I dr da 
Jo 



12(3 



+Y^v„,i?„,„3(r(r)rrr(o^))p(^) 



+ ^V„V;3K [dr [da 
^ Jo Jo 



,(0) 



v„,K dr (r ■■■e 



(2.42) 



where 0/2 is the Seeley-DeWitt coefficient computed with vanishing potential Vq. It 
can be read off directly from (p.38|) 



.(0) 



12 P 12 ""^24^ 1728 



uR. (2.43) 



One thus gets for the relevant terms we wish to consider 



^^^y<i V'^i?( Ji - J2 + J3) + ^R,u V^V'Vqi.h + J5 - 2 Je + 2 J7) 



+^RnVq MJu - J12 + 2Jio) - 



+^VqnR 




12 ^ 2 ^ 4 ^ 288 



+ ■ • • (2.44) 
where Jj are the integrals reported and calculated in appendix |A.3| . We finally obtain 



(2.45) 
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where 



6 V 12 2 4 ^ 288 / ^ 
1/^2 5„ 2„ 1 



is the "current" whose divergence gives the total derivative term. 

It is then clear that ( p.45| ) yields an unambiguous p-independent partition func- 
tion Z{(3). For p(r) = 6{t) one has the DBC propagator which yields directly 
the correct terms belonging to the local Seeley-DeWitt coefficient 03 contained in 
P) mi]- Others p's must then yield the same local result up to a p-dependent 



covariant total derivative^. Thus, we tested successfully the use of RNC to obtain 
the correct covariant results for Z{f3) also in the presence of external potentials. 

3. Supersymmetric nonlinear sigma models 

The N = 1 supersymmetric nonlinear sigma model is relevant for the worldline 
description of a spin 1/2 particle coupled to gravity. We have described this worldhne 
approach in [Q, where the representation of the one-loop effective action was written 
in terms of a path integral with periodic boundary conditions (PBC) for the bosonic 
coordinates x'^ and antiperiodic boundary conditions (ABC) for their supersymmetric 
partners, the worldline fermions ip'^. The worldline fermions play the role of the 
gamma matrices of the Dirac equation and the ABC are necessary to compute the 
trace in the spinorial space on which the gamma matrices act. These boundary 
conditions break worldline supersymmetry. In particular, the antiperiodic fermions 
ip'^ do not have any zero mode. The zero modes of the bosons can then be treated 
as in the previous section. 

On the other hand, PBC for both bosons and fermions preserves supersymmetry. 
The change of boundary conditions from ABC to PBC for the fermions corresponds 
to an insertion of the chiral matrix 7^ inside the trace in spinorial space. This gives 
directly the regulated expression of the Witten index for the nonlinear sigma model 
p!5[] , which is identified with the chiral U{1) anomaly corresponding to a massless 



Dirac fermion coupled to gravity |T^. In the present case, the periodic fermions 
acquire zero modes as well, and it is of interest to consider their factorization. In this 
section we describe this factorization using directly the action written in components. 
Superspace methods may also be used, but they do not seem to bring in drastic 
simplifications. 

^The coefficient 03 is related to the trace anomafies in 6 dimensions [E2|, which have been 
calculated also using the quantum mechanical path integral with DBC in p3[. 
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The = 1 supersjTiimetric nonlinear sigma model is described by the action 

^[X,^] = liy^g^,(x){x^X-' + r^) (3.1) 

where = ip'^ + x^T'^p{x)ipP is the worldline time derivative covariantized with the 
target space connection. We consider PBC for both bosons and fermions, x^(0) = 
x'^(l), ip'^i^) = 'ip'^il)- The corresponding partition function (Witten index) is then 



VxVij e-^[^'^] . (3.2) 



In the flat limit {g^i, = d^j^^) zero modes appear for both x^ and %p^, corresponding 
to the invariances 5x'^{t) = and 6iIj^{t) = 9^ for constant e'^ and 9^. In curved 
target space these two invariances are generically broken. Nevertheless we wish 
to factorize these "would be" zero modes to be able to carry out a perturbative 
evaluation and at the same time present the partition function as an integral over 
them. Thus we proceed as in section 2, and introduce extra dynamical gauge variables 
together with suitable gauge fixing conditions. The introduction of extra gauge 
variables can be obtained by considering the following identities 

S[x, = S[xq + y, = S[xq + ?/(xo, i), 

= S[xo + y{xo, O^i'ii'o + X)] = S[xo, ^o, ^, x] (3-3) 

where we first introduce the new gauge variable Xg, then change coordinates to 
RNC centered at Xq (this change to new coordinates and ip'^ is specified by the 
functions y^{xQ,C,'^) and ip^ {ip")) , and finally introduce a new gauge variable ip^ by 
setting ip^ = ipQ + x^- We end up with the action S[xo, ipo, ^, x] for the sigma model 
in RNC 

S[xo, ^0, X] = ^ ly \g,.{x,, + « + x^) ^^^°^+^'^^ ) (3.4) 



where the metric g^uixQ, ^) in RNC is given in (|2.19|) . This action contains the desired 



shift gauge symmetries encoded in the following BRST symmetry, suitably extended 
to two pairs of nonminimal auxiliary fields 



6x^ 


= ?7^A 






-Q''.{xo,ar))r]'^A 




= 


5x^{r 


) = 






= 


S-f^ = 









= «vr^A 


STTf, = 









= iPf,^ 


SPt, = 


. 





(3.5) 
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The nonhnear bosonic gauge symmetry acting on (xq,^'^) is just as in the previous 
section, while the fermionic symmetry acts linearly on the {ipi^, x^) fields. We gauge 
fix by choosing the gauge fermion 

^ = % [ dTp{T)e{r)+% [ drp{T)x^{T) (3.6) 
Jo Jo 

parametrized by the single function p(r) normalized to drpij) = 1. This gauge 
fixing is covariant under reparametrization of , and yields the following gauge fixed 
action 

•1 rl 



^[xo,^o,e,x]+«vr^ / dTp{T)e{T)-% / c/rp(r)Q^(xo,e(r))r7'^ 







1 



-tp^ / dTp{T)x''{T)-%r ■ (3.7) 

The ghosts 7^, 7^ can be trivially integrated out, while integration over the lagrange 
multipliers vr^ and produces the constraints 

dr p{T)e{r) = , [dr p{T)x^{r) = . (3.8) 

These constraints permit to invert the kinetic term to find the perturbative propa- 
gators for the fields and x^- Thus we end up with the path integral 



X exp ( - S[xo, ipo, ^, x] - Sfp[xo, ^, r], r]]) (3.9) 



where Spp is the same Faddeev- Popov action of the bosonic model written in ( |2.33| ). 



Note that since we use fermionic fields with curved indices {ip^, x'^) their nontrivial 
path integral measure is exactly compensated by the nontrivial measure of their 
bosonic supersymmetric partners (xg, ^'^). Therefore the complete measure is fiat 
and no extra ghosts are needed in this C3iSG, clS shown in IHl. 

The path integral can now be computed. Its value should not depend on (3 
as consequence of supersymmetry [|l5l. Thus a semiclassical calculation produces 



already the complete result JT^. However, we have now set up the path integral in 
such a way that higher loop calculations can be unambiguously performed to test 
the (3 independence. We perform this test in the next subsection. 

3.1 Order P correction to the chiral anomaly 

In this section we consider the perturbative expansion for the = 1 supersymmetric 
sigma model with periodic boundary conditions on all fields, as described in the 
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previous section. We will explicitly verify the {3 and p independence of the path 
integral for the chiral anomaly Iw in D = 2 and D = 4 up to order /3 (i.e. two 
loops on the worldline), employing dimensional regularization |^ whenever necessary. 
Before starting, let us recall that the chiral anomaly Iw can also be written as 

Iw = Tr ((-1)^ e-'^^^) = Tr (75 ei^ (3.10) 

where Q = -^j^ is the supersymmetry charge (the Dirac operator), H = = 
— = — |-R the quantum hamiltonian, and (—1)'^ = 75 the fermion number 
operator of the N = 1 supersymmetric sigma model. Dimensional regularization 
preserves supersymmetry without the need of any counterterm 0: the quantum 
hamiltonian H must have a potential obtained from squaring the supercharge 
Q, and DR produces exactly that potential. Thus the path integral representation 
of Iw in (|3.9| ) is good as it stands when using dimensional regularization. Other 
regularization schemes may need counterterms to enforce supersymmetry. 

Let us consider the perturbative expansion at fixed Xq and ip^. One obtains the 
various propagators from the free action 

82 = ^ f^dr + X'r) - V.V' (3.11) 



with the fields constrained by ( |3.8| ). The propagators for and the ghosts r]^^r]^ 
were already described in ( p.26| ), whereas the one for the fermions is given by 

{x'{r)x''{<y))=Pg'"'{x,):F{r,a) (3.12) 

where the Green function 

J^(r,a) = •A(r - a) - f da' p{(j') •A(r - a') + [ da' p{a') *A{a - a') (3.13) 

Jq Jo 

satisfies 

'T{T,a) = 5{T-a)-p{T) 

-r{r,a) = 5{T-a)-p{a) (3.14) 

(dots on the left/right denote derivatives with respect to the first/second variable) 
and 

1 

c/r p(r) J^(r, a) = . (3.15) 







With these propagators the chiral anomaly is given by 

(27ri)T 



iw 
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where S'^*"*^ = Sg — Sq^2, Sg = S + Spp is the total action in ( |3.9|) , and (1) = 1 is the 
normahzation of the (p-dependent) path integraL For the sake of compactness, it is 
useful to rescale ipo — >• ^/P ipo, getting d^%jjo — >• P~^^'^ d^ipo and ip a/5 ipo+X- Being 
ip of order ^/J3 we can organize the interaction terms in a systematic way, namely 
we split Sq"'^^ = Sq^3 + Sq^4 + . . . wlth the Sq^n = Sn + Spp^n coutributious of order 
0{P^~^). Moreover, in the RNC expansion 5*3 = Sfp,3 = 0, as g^u,x{xo) = 0. Thus, 
after some simplifying manipulations that exploit index symmetries, the relevant 
vertices for our purposes read^ 

^4 = ^Rc.,up fdr + + xnr) 

+ ^i?a;..A f^dr C^^i^/Pro + xn{^/P^o + X') (3.17) 
^5 = ^5,1 + ^V^i?„^.A j^dr C^^e{^/Pro + xni^/Pi'o + x') (3.18) 

Sq = 5*6,1 + — (—Rxa/SaR'T^iu'^ + '^RxauaR/S^j.-y'^ + -^RxafiaR^u^T 



+ ^S/^VpR^^,^ jdr re^ee^V^^ 



P% + xnWP% + Xl (3.19) 

where all tensors are evaluated at the point xq. 

It is clear from the Grassmannian nature of the measure that the only terms 
that can give an a priori nonvanishing contribution to ( p.l6| ) are those that saturate 



the measure itself. Hence, in D dimensions only terms with D fermionic zero modes 
can contribute; in particular 



d''i)o ■ ■ ■ = e^^-^"" . (3.20) 

Before starting our calculation, let us briefly recall how some known results for chiral 
anomalies in arbitrary dimensions D can be derived from ( |3.16| ). First of all, the 



chiral anomaly is trivially zero when D = 2/c — 1, G N, because an odd number 
of ipo is needed to saturate the measure whereas the fermions always come in pairs. 
Conversely, if one considers D = 2k, it is not difficult to convince oneself |T^, 



that, at the one- loop level, the only nonvanishing contributions come from the k- 
th power of the vertex \Ra^j.uX Jq dr ^""^^ipQipQ ; hence, no fermionic propagators are 
involved. More specifically, the independent terms contributing to the anomalies are 
of the form 

Tr [n'] Lk , n^p = Ra(3f.u <V^o (3.21) 



'^s,! and 56,1 are contributions to ^5 and 5*6 that come from the cubic and quartic order in the 



expansion of g^.v{x) in RNC, see (2.19). They are not exphcitly reported here since they will not 



enter in the forthcoming calculations. 
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where the integrals Lk are given by 



Lk= [dn--- [drk'B{n,T2)'BiT2,T,)---'BiTk,n) . 

Jo Jo 



Integrating by parts and making use of the identity *B* 
reduced to their string-inspired counterparts 



(3.22) 

*A*, the latter can be 



dri 



dn •A(ri - T2)'A{t2 - n) ■ • ■•A(rfc - n] 







C(2/) k 



(27rj)2' 

k 



21 

21 + 1 



I e N 



(3.23) 



where ( is the Riemann zeta function. Notice that for k odd, both traces and integrals 
vanish: indeed, gravitational contributions to the chiral U{1) anomaly are present 
only for spacetime dimensions D = 4n, G N. The result of this semiclassical 
approximation is known to be exact |]14|, and it is independent of the boundary 
conditions chosen for the quantum fluctuations (see also p4|). Nevertheless, the full 
path integral ( p.9| ) should be independent of P and p. We have explicitly tested this 
property by doing a two-loop calculation for the simplest cases D = 2,4. For D = 2 
each worldline graph vanishes, so the proof is quite trivial. Thus, we directly pass to 
describe the = 4 case. 

The two-loop truncation of ( |3.1(j| ) reads 



"(WP 



exp 



(3.24) 



where we have omitted those terms that are explicitly zero; the subscript c refers to 
connected diagrams. To further simplify this expression, let us make few preliminary 
considerations. Recalling the splitting Sq = S + Spp, we compute 



(54) 



^Ra,up fdr (re^(er + (V^< + 

6p Jo 



dr {cewpro + xnwPi^o + X')) 



(3.25) 



where we have used some of the formulae in appendix |A.2| . This term does not 
contain fermionic zero modes. This is also true for all the vertices coming from the 
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ghost action Spp- We also note that (5*6) is of order (3^ and contains at most two 
zero modes. This imphes that 5*6 can only enter at two-loops through connected 
terms. Thus, the (one-loop) anomaly is given by 



r(0) 



(27r)2 2 ■ 42 



(3.26) 



where from eq. (|3.23| ) L2 = — ^, as expected. 

From these preliminaries one recognizes that the order (3 correction is given by 



+-^{Sfp,a){S'1)c + -{Sfp,4,SI), 



(3.27) 



We first compute the four terms in the first line, leaving the contributions from the 
FP vertices in the second line for the very end. 

The first term can be immediately obtained from (|3.25|) and ( |3.26D . Using the 
identities listed in appendix |A.2| , it can be written as 



1 f d^XQ d^ipo 

(271/3)2 



{S4){Sl)c 



(3 f d'^XQ 



384 y (27r)2 



C, i2Ko - SK2) . 



(3.28) 



As for the second term, the only way to pick out four zero modes is to square the 
connection part, the one explicitly reported in ( p.l8| ), obtaining 



2 \ (4^-0) 



;si> 



36 



X [ dn [ dT2 (C'^'^'i^' ^"2^/32^M2 

Jo Jo 



36 

+ V„i?^<,AiA2(V^/?"^3A4 + V"i?'^^3Aj(/5-/6)] • (3.29) 



Using the identities in appendix A. 2 it becomes 



1 f d^xo d'^ipo 



[SI), 



d'^Xo 



--{Ko + 4K,) (2/2-/3-/4) 



+ ^K2 (4/2 - 2/3 - 2/4 + 3/5 - 3/0) + ih - h) 
18 6 



. (3.30) 
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The third term yields 



(4Vo) 



Jo Jo 



.^Ml6(^P2-2^/3l72 ^ ^P272^/3l-2 ^ ^0272^/3lP2) (3 3^) 



where the six-index tensor is the one given in the round parenthesis of equation (|3.19| ) 
After a bit of tensor algebra one gets 



(2vr/3)2 

Finally, we are left with the fourth and more involved term 



. , —Ko + -Ki - -K2 + -A^3 ) h ■ (3.32) 
(27r)2 V 24 " 9 9 3 ' ^ ^ 



(5|>f^°) = ((S4,l + ^4,2f> 



(3.33) 



where 5*4,1 and 5*4,2 can be read off from ( |3.17D as the first and second term, re- 
spectively. The former may yield only one zero mode and thus it does not give any 
contribution to ( |3.33|) by itself. The remaining terms give instead 



3 f d xq d ipo 



3! J {2TTpy 
3 f d'^xo d^ipQ 



3! j (27r/5)2 



((5*4,1)^5*4,2)^ 
(5*4,1(5*4,2)^)^ 



(3 f d'^xo 



144 J (27r)2 



-/3 



d'^Xo 
(2^ 



-Ko (3/i4 + Is) 



(3.34) 
(3.35) 



+—K, (-2/11 + 4/13 - 2/12 - Jio - I9) 



+ ^K2 (6J11 - 12/13 + 6/12 - 6/8 + 3/10 + 3/9) 



1 f d^Xp d'^iPo , ^3\ _ P f d^o J J 

3! j (2vr/3)2 - 16 y (2^0^ ~ ~ ^''^ 



(3.36) 



We now explicitly use the integrals listed in appendix [A.3| , so that the four terms 
computed so far sum up to 



d^Xp 



288 



12 



72 V 36 3^ 



1 



+ 72 I'S^'^ + ^^^^liij^' + M^^"^^^' ' 720 



720 
13 



(3.37) 
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This is non zero unless one adopts DBC Feynman rules, for which Cdbc = ~j2^ 

1 f ^ 

Ebbc = tIo' so that I^^ oc / ^^Ki = by Bianchi identities. However, the 

1=1 

Faddeev-Popov action gives the additional contributions 

1 f d^xo d^^ipo 



(3 



d^XQ 
(2^ 



[(5^P,4)((54,2)'>,+ (^FP,4(V)'>J 



(3.38) 



which lead to the complete order-/? correction 



-P 



d'^Xo 1 
(27r)2 72 V36 



13 
720 



3 

E 

1=1 



Kj = . 



(3.39) 



This result shows that one may safely adopt p-dependent Feynman rules also in the 
case of periodic worldline fermions, provided that one appropriately takes care of the 
zero modes. As we showed, this can be done through the appropriate introduction of 
auxiliary gauge variables, together with a suitable gauge fixing of the ensuing shift 
gauge symmetry, which turns the auxiliary gauge variables into the wanted zero 
modes. As a consequence of this procedure, one must include the corresponding FP 
determinant, which vanishes only in the DBC case. Indeed, in where a similar 
calculation has been carried out in DBC as a test for the time slicing regularization, 
no extra vertices were needed to obtain the correct result. 



4. Conclusions 

We have discussed in great details the factorization of perturbative zero modes which 
appear in the worldline approach to one-loop quantities for spin and spin 1/2 
particles. We have focused particularly on the issue of reparametrization invariance 
of the corresponding nonlinear sigma models. As well-known, the calculation of 
one-loop quantities in the worldline formalism requires to master the path integral 
for sigma models on the circle. In perturbative computations around flat target 
space, perturbative zero modes appear. We have described their factorization using 
BRST methods and employed an arbitrary function p{t), the so-called background 
charge, to parametrize different gauges. The arbitrary function p(r) allows to test 
for gauge independence and includes as particular cases the two most used methods: 
i) The method based on Dirichlet boundary conditions (DBC), which is directly 
related to the calculation of the heat kernel. It produces covariant results in arbitrary 
coordinates also before integrating over zero modes, but it has a propagator that is 
not translational invariant on the worldline. ii) The "string inspired" method, where 
the constant modes in the Fourier expansion of the quantum fields are factorized. 



- 23 - 



It yields a translational invariant propagator, but when using arbitrary coordinates 
produces noncovariant total derivatives in the partition function density. The string 
inspired propagator is simpler and faster for computations on the circle, and the 
technical problem of the noncovariant total derivatives can be overcome by using 
Riemann normal coordinates. 

The BRST symmetry guarantees the gauge independence of the partition func- 
tion Z{(3). From this one deduces that the partition function density z^^\xo, (3) 
should contain a universal p-independent covariant part and at most a p-dependent 
total derivative, which vanishes upon space-time integration. The p-dependent total 
derivative is in general a noncovariant expression when calculated in arbitrary co- 
ordinates, as exemplified in eq. ( p.lO| ). However, it is covariant if calculated using 
Riemann normal coordinates. The source of this discrepancy is easy to understand: 
the gauge-fixing constraint is expressed in terms of coordinate differences 
and does not transform covariantly as a vector under reparametrization of the back- 
ground point Xq, but it does once one chooses Riemann coordinates centered at Xq. 
The gauge fixing in Riemann coordinates is rather subtle. It has originally been 
studied by Friedan for applications to bosonic 2D nonlinear sigma models 0, and 



recently employed by Kleinert and Chervyakov |T0[ in the ID case. It makes use of 
a suitable nonlinear shift symmetry. 

In this paper we have first addressed some issues left open in and performed 
extensive tests to check the correctness of the BRST method for extracting the zero 
modes. Then we have extended this method to the supersymmetric nonlinear sigma 
model, which is relevant to the description of spin 1/2 particles coupled to gravity. 
All our tests have been successful. We can now draw the following conclusions: 

• The perturbative calculation oi Z{P) is gauge independent if the external fields 
describing the interactions have derivatives vanishing sufficiently fast at infinity. This 
is a necessary condition for the validity of the perturbative expansions employed in 
our tests. If these conditions are met, all total derivatives integrate to zero. When 
the external fields have derivatives which do not vanish fast enough, one may still 
proceed for example by trying to resum a class of diagrams and make the perturbative 
expansion of Z{(3) well-defined. This works in the harmonic oscillator, which has 
the potential V{x) = ^u'^x'^. In this case we have explicitly checked the gauge 
independence of the partition function up to order (Puj)'^. 

• The covariant local expansions of z^^\xo, P) in Riemann normal coordinates 
(RNC) are different for different p's. The difference is given by covariant total deriva- 
tives with coefficients depending on p. 

• One can use RNC with any p to compute local trace anomalies in D = 4. In 
particular, one can use the string inspired propagator. The total derivatives only 
affect the nR term, which anyway is a trivial anomaly (it can be canceled by a 
counterterm) . The reason why the "string inspired" computation tested in failed 
is because the Faddeev-Popov contribution arising form the nonlinear shift symmetry 
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needed for the implementation of RNC was missing. 

• The use of RNC yields local covariant results also in the presence of external 
fields like scalar and vector potentials. We have explicitly verified this in the presence 
of a scalar potential V. We have performed a higher loop calculation to make sure 
that the expected results are indeed obtained. It is not obvious how to use arbitrary 
coordinates to check this property. In particular, the method proposed in section 7 



of [10 1 to achieve a local covariant results in arbitrary coordinates does not seem to 
be correct, as for example it does not yield the correct covariantization of the higher 
derivatives of the scalar potential. 

• We have used the supersymmetric model with the correct factorization of the 
zero modes to test successfully for any p the P independence of the Witten index, 
i.e. of the chiral anomaly for a spin 1/2 particle. We have computed the order (3 
corrections the chiral anomaly in D = 2 and D = 4 and found that they vanish. 

The results described here make sure that path integrals for nonlinear sigma 
model on the circle, with and without supersymmetry, are in good shape both in 
their UV and IR structure. They can be employed to produce unambiguous results 
in the worldline formalism with background gravity. The method discussed here 
could also turn out to be relevant for the extension to the noncommutative case, 
recently studied in |^ for describing the coupling of -DO branes to gravity, and to 



the free field analysis of the AdS/CFT correspondence [p7 
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A. Appendix 

A.l Linear sigma model 

Let us consider the linear sigma model 

'1 



1 



(A.l) 



and compute perturbatively in jd the partition function 

Z = Ti e"^^ = i Dx e-^[^] . (A.2) 
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We first try to understand the factorization of the zero mode with some generahty. 
To extract the zero modes one can first introduce a shift symmetry, and then gauge 
fix it in a suitable way. Using the Faddeev-Popov method one obtains the following 
chain of identities 

Z ^ j Dx e-^[^l 

Dy e-^I^o+^'l 
^ ^ Dy J d^e J dTp{T)y^{T)^ e'^l^o+y] 



S[xo€+y] 

= j d'^xo fDy5(^j drp(r)y^{r)^ e-^i^^o+y] ^ J ^(p) (^r,^^ . (a.3) 

Here we have first used the linear split x^{t) — Xq + y^ir) for an arbitrary constant 
Xq , and then the translation invariance of the path integral measure Dx = Dy. Then 
we made use of the shift invariance (typical of the background field method) 

•*'0 -^Oe — Xq -r t 

y^^ ^y^ = y>^-e^ (A.4) 

which leaves invariant the field x^{t) since x'^(r) = Xq +y^{T) = x^^ There- 
fore the action S[x] itself remains invariant. One may then use the Faddeev-Popov 
trick of inserting unity represented as follows 

1 = J d''eS^drp{T)yi:{T) 

= J d^'eSi^j dTp{T){y'^{T)-e^)) 

= J d^eS(^J dTp{T)yi'{T) - e") (A.5) 

where the background charge is normalized to unity, dr pir) = 1. Finally we used 
that d^e = d^XQ^. This formally proves that the final result for the partition function 
Z = J d^XQ (27r/3)~T^(/')(a;o) is gauge invariant, i.e. independent of the choice of the 
function p. However, the density z^'^\xo) can in general depend on p through total 
derivative terms which should vanish upon integration over the zero modes Xq. 

One can derive the same result using BRST methods. The action S[xo + y] can 
be considered as a functional of both Xq and y^{T) and is invariant under the shift 
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]4D. To fix this gauge invariance one can introduce a constant ghost field r]^ for this 
abehan shift symmetry. The corresponding BRST symmetry transformation rules 
are 

Sy^" = -ry^A 

5r/^ = . (A.6) 

To gauge fix one must also introduce constant nonminimal fields fj, it with the BRST 
rules 

5tt^ = 0. (A.7) 
m=r^JdT pW{t) (A.8) 

^0 



Using the gauge fixing fermion 



we obtain the gauge fixed action 



5 



Sgf = S[x^ + y] + —^ 



S[xo + y]+ iTCf, / dr 
Jo 



p{T)y'{r)-n,v' (A.9) 



where ^ denotes a BRST variation with the parameter A removed form the left. The 
ghosts can be trivially integrated out while the integration over the auxiliary variable 
TT^ produces a delta function. Thus, the BRST method reproduces the last line of 
( [A.3|) , however it makes it easier to implement more general gauges (for example to 
get massive propagators), though they will not be considered here. 

To test the previous set-up, let us compute perturbatively in the proper time (3 
the partition function in eq. ( |A.2|) . From (|A.3| ) one can write it as 



D_ 

2 



z^p\xq) = (e"^'"*) = exp[(e-^""')c - 1] (A.IO) 

where (27r/5)~"? is the well-known normalization of the path integral measure, and 
the quantum average (e~'^™*) is obtained using the free propagators (i.e. from the 
action with = 0) of the fields constrained by the gauge fixing. These quantum 
fields are given by y^(r) = x^{t) —Xq with Xq = dTp{T)x^{T) and their propagator 
reads 

(y"(^)y^W) = -/55'^^S(p)(r,a) (A.ll) 
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where the Green function i3(p) (r, cr) has already been described in eq. ( p.28 



We now aim to compute eq. (|A.2|) to order by expanding the interacting 
action around the constant Xq 
•1 



S^nt = 1^ dr /?V(x(r)) 
P Jo 



= ^ fdr (^ + \y^y''dpd,V + ■ ■ ■ + ^l/'^V +...) (A.12) 

"^^ •S'e S\o 

where for Siq we used an obvious short-hand notation. Here all vertices are evaluated 
at Xq (i.e. one sets x^[t) = Xq + y'^ij) and expand the action around the constant 
Xq). The vertices indicated by S'„ contribute effectively like t/". Considering that 
5*4 gives only disconnected contributions and that only an even number of y give 
nonvanishing Wick contractions we obtain 



= exp 



— (5*4 + 5*6 + 5*8 + 5*10) + 2^^'^ + 'S'e + 25'5S'7)c 



(A.13) 



and thus 



Z = ^ / rf^xo exp \-/3V+ ^M,dV - ^ f^Marf^ + M,{d^Vy) 

(27r/5) 2 J VI 2 V4 / 



+0(/?' 



(A.14) 



where the Wick contractions obtained using ( [A.ll| ) produce the following integrals 



Ml 
M2 
Ms 
M4 
Ms 
Me 



rfr i3( 







1 

dr i32(r,r) 

1 

dT B%T,t) 

1 .1 

dr / (icr i3(r, cr) = Cp 

Jo 

1 .1 

c/r / (icr i3^(r, cr) 



1 dr da B{t, a)B{a, a) 
Jo 



1 

12 



Co + 2c: 
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where on top of Cp given in eq. (PI) we have defined C; = drF^ir) and C'^ = 
/o drF^ir). Recall that the DBC propagator is obtained by setting p(r) = 6{t) and 
the SI one by p(r) = 1. In the SI case all these Cp's vanish, while the particular 
values of the integrals for the DBC propagators are as follows 



Ml = , M2 = 4t , ^3 



1 1 
M4 = , M5 

140 ' 12 ' ^ 



90 ' 60 



Let us now consider the partition function at order (3^ (we rename xq by x) 
1 



d^x 



I - f3V + ^{V^ + MiuV) 



+^l-(-4V^3 _ i2MiV^nl^ - SMsdV - 12M4(a^\/)2) 



+0(/3^ 



(A.15) 



From this we deduce that 

AZ = ^(arbitrary p) - Z{DBC) = ^ [ d^x\ — {Mi + -)nV 

(27r/?)T J 12 6 

/?3 



24 



3(M2 - 4)°'^ + 12(Mi + i)9^(V5^y) 
-30 o 



+1{M, - ^)d^{d''Vd,d^V) + i(M2 - ^)d^{Vd,nV) 

+hM, + ]:)d^{V%V)\ +0{f3') 
4 J 



(A.16) 



is a total derivative. For potentials V with derivatives vanishing sufficiently fast at 
infinity (a condition which is necessary for the validity of the perturbation expansion 
in P) the total derivative can be dropped and all ways of factoring out the zero modes 
are equivalent. 



A. 2 Conventions and identities 

We use the following conventions for the curvatures 



(A.17) 
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The change of coordinates to Riemann normal coordinates is given by 

y'{xo, = e~Yl -^"3 ■ ■ ■ v„„r^,,,(xo) r--- (A.18) 



n=2 



where the covariant derivatives act on lower indices only. 

We have found it convenient to express all the terms in the order-/? correction 
to the chiral anomaly (section |3.1|) as combinations of the invariants 

T/- 01---04 po- p TD P^"^ 

Ty^ pCn---04 po" p p f?^"^ 

J?- _ p«l - «4 po- p p TD/liy 
-'^S ^ -'I' ^ai ^^crvpa2 0304 5 

where i^i + K2 + -ft's = because of Bianchi identities. By using Bianchi identities 
and integration by parts one can easily prove the following identities 



j d^xe''^-'''RR^,o.,a,R^\,a, = j d^x {2K0-8K2) 

which have been used to cast the final results in a more compact form. 

We now list few useful identities involving the propagators. Recalling the defi- 
nition of A(r — cr) in ( p.29| ), we start by defining the quantities 



Cp= dn / dT2 p{ti) p{t2) A(ri - T2) [Cdbc = -j^] 
°i % 1 

Dp= dn / dT2 / drs p(ri) p(r2) p(r3) A(ri - T2)A{ti - t^) [Dbbc = ni] 
Jo_^ Jo^ Jo 

" lo^^^Io^^^ ^^^^"^ ^^^^^ ~ ^^'^ ~ 7^ " 

which all vanish in the string inspired approach (p('r) = 1). For the p-dependent 
propagators we need 

B{t, a) = A(r - a) - [ da' p{a') A(r -a') - [ da' p{a') A(a - a') + Cp 

Jo Jo 
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B\, = B{r, r) = C, - ^ - 2 J^da'pia') A(t - a') 

•B{t, a) = •A(t -a) - [ da'p{a') •A(t - a') 

Jo 

•B\r = V3{t, t) = - fda'p{a') •A(t - a') 
Jo 



-^B\r = 2'B^r 



dT 

•B'{T,a) = •A*(r-a) 

J^(t, a) = •A(r - a) - [ da'p{a') 'A{t - a') + [ da' p{a') •A((j - a') 

Jo Jo 

T\-r = J^{T,r) = 

-B'{T,a) + 'J'{T,a) = l-p{T) 



and the identities 



1 

dw •A(t - u;) •A((7 - a;) = -A(t - cr) 



These identities have been used to the express aU the worldhne integrals, hsted in 
the next appendix, in an economical form, namely as combinations of Cp, Dp, Ep and 
pure numbers. 

A. 3 Integrals 

We list here the worldline integrals needed in the main text, evaluated with di- 
mensional regularization whenever necessary. We use the notations B = B{ti,T2), 
B\t = B{t,t), and similarly for its derivatives. 

In section 2.1 (partition function at 3 loops), we needed the following integrals 



[drB\r CB' 
Jo 



IdrBl (•B* + A,,)|, = C, 



1 

12 



1 

H,^ ldTV3\l^-Cp 





Hs = p{t) B\r^-Cp-^ 



H,= IdT B\l i'B' + A,,)|. = - -Cp -^Ep + — 
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1 /•! 

2 



i/g - / / dT2 -fill -B- B-|2 B = - -C, - 



^0 

1 rl 



H, = jdr, jdT2 -^li 'B B-\2 B' = Cl + ^C, -D,+ ^E, 
Hio = J^dnjy2 'B\, 'B- B\2 B- = Cl - ic, -D^- ^E, 
Hu = J/^^jy^ *^ ^ (*^* + \h)\2 = -Cl + + + ^-E, 
Hi2 = J}^i f}^2 B\, 'B' {'B' + A,,)|2 = -Cl + ^C, + + 1^^ - 

^13 = j^^j}^^ C^" - ^Ih) = -Cl + + + 



11 

1440 



Hm = [dn [dT2 'B' B'^ = Cl -2D, + Ep + ^ 
Jo Jo 

Hi5 = 1^^!^'^ ^' cl - ^Cp -Dp- ^Ep - 

//i6 = fdr p{t) B\l = -?>Cl + ^C, + 4L>, + ^ 
H,, = J}^^jy^ p(ti)p(t2) = - 2L>, + + ^ 

Hi, = l'dnJ^dT2 p(t2) 'B' B\i = - Ac, _ _ _ -i 
i^2o = J/^^jy^ S -fill = cl - ^Cp - ^Dp + ^E, 



720 

144 



In the 4 loop calculation of section 2.2 we made use of the following integrals 

Jl = jy^j}^^ ^ C^* + ^9h)\2 = Cl - ^Cp - ^Ep 

J2 = J}^iJ/^^ ^ = -Cl + yCp + \dp + 1^;, 
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J3= [ dn [ dT2 p(r2) B B\2 = -20^ + 2D, 
Jo Jo 



720 
1 

144 



720 

P 18"'' 3^^ ' 144 ■ 



J^ = j\TB\ = Cl-^C,-lE,+ ^ 
Jg = j'dT B\^ = Cp-^ 



1 

Jio = I dr p{t) B\r = -Cp - ^ 



Jn = ^rfr B\^ {'B' + Ag,)|. = ^^'1^ 

Jl2= CdT'B\ = -Cp. 







Finally, the integrals needed in the calculation of the order-/? correction to the chiral 
anomaly in D = 4 (section |3.1| ) are given by 



h= j dnj dT,_B\rB-B\-i = -2h 

Jo Jo 

h = [dn fdT2 -Sli B 'B\2 = ^Cp + \Ep 

Jo Jo 36 6 

h = [ dn [ dT2 BU'B' B\2 = 

Jo Jo 

h = ! dn f dT2 *BB* B = -\Cp -\Ep+ ^ 



Jo 



^ ' Q ' 360 

>i /.I 

/g = I dn dT2 B B'B' = -2/5 



Jo 



7 



/ dr, [ dT2 [drs {B •^•)(i,2) i3(2,3) '^(3,1) = h 
Jo Jo Jo 
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'0 



dri 



dri 



[ dT2 f drs B\r:F^i,2) (•BB*)(2,3) = 

Jo Jo 

dr2 [ drs {B •^)(i,2) '^(2,3) •%,!) -- 

Jo 

dT2 I dn {B* *T\^,2) B(2,3) 'Bo,!) : 



^0 



Cn — ~En 

9 " 3 

-Ir - If - J- 

36 " e, " 720 



^0 



^0 



JO 



^0 



'dT2[dTs{'B' + 'J')UBl,^ %3) 

Jo 

\t, [dn BU -^(1,2) -^(2,3) ^(3,1) = 
dr2 fdn 'BU B(l,2) -^3^2,3) ^(3,1) = 

Jo 

\t2 CdT^ {*B B*\,,2) {'BB-\2,3)-- 

Jo 

'drJdr.iJ^ i3-)(i,2) %3) ^(3,1) = 
Jo 

dT2 / rfrs (J^ B)(i,2) 'iS'a^a) ^'g^^) = 

Jo 

fdT2 fdn (-F •i3-)(l,2) i3^2,3) ^(3,1) 

Jo Jo 

[ dn [ dT2 [dT3 p(Ti) S(i_2) -6(2,3) -^CS,!) ^ 

'o Jo Jo 



'■6(3,1) = —Co + T^Eo 



18 



Jo 



Jo 



12 ' 


1 

+ 

144 










36 


6 ^ 


1 

^ 144 


36^" 


-\ 

6 


1 

^ 360 


36 


3 " 


1 

^ 360 


18 ^ 




1 

' 720 


^^^"^ 




1 

720 ■ 
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